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ABSTRACT 

The velocity field in the lower solar atmosphere undergoes strong interactions with magnetic fields. Many authors have pointed out 
that power is reduced by a factor between two and three within magnetic regions, depending on frequency, depth, the radius and the 
magnetic strength of the flux tube. Many mechanisms have been proposed to explain the observations. In this work, SDO dopplergrams 
and magnetograms of 12 bipolar active regions (jSARs) at a 45 second cadence, are used to investigate the relation between velocity 
fluctuations and magnetic fields. We show that there is an asymmetry within jSARs, with the velocity oscillation amplitude being more 
suppressed in the leading polarities compared to the trailing polarities. Also, the strongest magnetic fields do not completely suppress 
the five-minute oscillation amplitude, neither in the spot innermost umbrae. 
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1. Introduction 

The study of the interaction between waves and magnetic fields 
plays an important role in the investigation of the mechanisms 
responsible for the solar corona heating. 

It is well established that the amplitude of solar five-minute 
oscillations in the phot osphere is reduced by a factor of two 



to three within 



1981 


; Lites et al. 


1987 


; 'Tarbell et al.' 


tHindman & Brown 



ARs (lLeightonetal.Ml962 



Abdelatif et al 



-r . Woods & Cram 



Braun et af. 
Title etal.|ll9'92[ 



1988; ' Braun et al.lli988F 
1998). The reduction is observed to 
depend on frequency, depth, and on the physical parame- 
ters that descri be the flux t ubes, i .e. their radius and mag- 
netic strengfli jBrown et al.l 119921; iHindman & Brownl 119981: 



iGordovskvv & Jai'nr2008 ). This behavior has bee n also found 

concentrations 



in small 
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1996; 
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Hindman & JainI 2008; iJain et al. 


2009; 
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It is also known that at higher frequencies, above the acoustic 
cutoff frequency in the low photosphere, the velo city oscilla- 
tions a mplitude is enhanc e d at the edges of ARs (iBraun et al.l 
[T99I iBrown et al.1 IT991 iToner & Laboiitel [T99l . Also, an 
enhancement of power in the three-minute ba nd has recently 
been observed in the inner umbra of a pore (IStangalini et al.l 
I20T2I) . 

Sever al mechanisms have b een proposed to explain the observa- 
tions (IHindman et al.lll997h : 

1) The intrinsic power inhibition due to local con- 
vection suppression. B y usi ng numerical simulations 
iParchevskv & KosovichevI (l2007t) found that the reduction 
of wave excitation in a sunspot, can account up to 50% of the 
power deficit. 

2) Pa r tial p-mode a bsorption (iBraun et al.|[r987t Bogdan et al 



[T99I ICailvl IT995I ICallv & AndriesI 120101) . ISpruit & Bogdan 



interpreted in terms of p-modes mode conver sion between the 
oscill ations and MHD waves. Both models (' Callv & BogdanI 
1 199^ and simulations (Cameron et al. 2008) have shown that 
mode conversion is able to remove a significant amount of 
energy from the incident helioseismic wave, with an efficiency 
that depends on the angle of the magnet ic field from the 
vertical, with a max imum at about 30 ° ( ICallv et al.1 120031; 
,Crouch & Callvll2003l) . 

3) Opacity effects. Within a magnetic region, the line-of-sight 
optical depth experiences a depression (Wilson depression). Due 
to increasing density, the amplitude of the velocity fluctuations 
decrease with depth. 

4) Alteration o f the p-mode eigenfunctions at the ha nds of the 
magnetic field dJain et al.ll 19961 Hindman et alJll997l) . 

We note that the majority of the studies on the interaction 
between solar oscillations and magnetic field have focused on 
the use of models and simulations. 

Because of the dependence of the velocity amplitude on 
magnetic strength, the most noticeable effects are observed 
within ARs, where the strongest magnetic fields are found. 
Among them, bipolar active regions (hereafter ySARs) are 
particularly interesting, showing asymmetries in their mor- 
phology and physical properties, as s everal works have 
already pointed out (se e , e.g., Brav & Loughhead 11979 ; 

1980; Temulloetal. 1981; Zwaan 1985; 

Peti-oyav et al., ,1990 ; 
1994 conclude that 



Balthasar & Woehl 
van Driel-Gesztelvi & Petrovav, 



■1990^. 



Fan et al. 1993 ). Moreno-Insertis et al.l 



(I1992I) showed that the wave absorption by sunspots can be 



morphological asymmetries are due to the different inclination 
with which the pol arities of jSARs e merge. 
Using simulations. iFan eFal 1 (119931) have studied area asymme- 
tries, and predict that the leading polarity must have a strength 
about twice larger than the trailing polarity, at the same depth. 
This result, which is a consequence of the Coriolis force action 
on the magnetic structure, is able to account for the trailing spot 
fragmentation, its lower flux magnitude and shorter lifetime. 
MDI observations of 138 bipolar magnetic regions have recently 
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quantitatively shown that the areas of leading po larities are 
typica lly smaller than those of trailing polarities (lYamamotd 
I2OI2I) . This area asymmetry could be produced by the Coriolis 
force during a flux tube's rising motion in the solar convection 
zone. 

In this work we investigated the velocity oscillation ampli- 
tude reduction by the magnetic field within jSARs. We analyzed 
12 /3ARs from SDO-HMI data, both in the Nothern (N) and 
Southern (S) hemispheres, and we found that the leading po- 
larity systematically has a greater reduction in amplitude than 
the trailing. In addition, we studied the amplitude reduction as a 
function of the field inclination for ARl 11 66. 

2. SDO-HMI data 

Our data set consists of 12 SDO-HMI magnetogr ams and dopp- 
lergrams pairs with 1 ar c sec spatial resolution (IScherrer et al.l 
I2OI2L ISchou etani20T2h . These were acquired with a 45 sec- 
ond cadence, which set the Nyquist frequency. Each /3AR data 
set is 3 hours long, which set the lower cut-off at ~ lO"'*^"'. 
Observation times cover the interval from 201 1 March 8 to 201 2 
January 3 (see Table 1). 

We selected 8 isolated bipolar regions in the N hemisphere, and 
4 in the S hemisphere. In order to limit the effects due to the in- 
clination with respect to the line-of-sight, we selected /3ARs as 
close as possible to the disk center The rectangular area enclos- 
ing the ySARs was selected such that it included all the magnetic 
features with strength above 500 G, with the the minimum and 
maximum allowed dimensions of the enclosing rectangle being 
500 and 700 arcseconds, respectively. Figure [1] shows the mean 
magnetograms of two selected ySARs: AR 1 1 166 in the N hemi- 
sphere, AR 1 1316 in the S hemisphere. In both panels the lead- 
ing polarity is on the right and it is negative for the Northern 
ySARs, positive for the Southern jSARs. Table 1 shows the list of 
ySARs analyzed in this work. 

3. Method and analysis 

The selected regions around each /3AR were co-registered us- 
ing a FFT technique ensuring sub-pixel accuracy. Inspection of 
the co-registered data shows that the strong magnetic structures 
remain in the same position during the entire duration of the ac- 
quisition. The amplitude of velocity oscillations was estimated 
pixel-by-pixel through their integrated spectrum. Specifically, 
we selected a spectral window Av centered at 3 mHz (five- 
minute band), having a width of 1 mHz, and then we integrated 
the FFT amplitude spectrum Ay in that frequency range 



Jav 



Aydv 



We then divided the data into bins of magnetic field, each 25 
G wide, and averaged the amplitude of velocity oscillations in 
each bin, thus obtaining =< Aw >b- We consider only pixels 
with |B| > 25 G. Since all the amplitude distributions within the 
respective magnetic bin have a Gaussian-like shape (see inset in 
FigurelSJ)), we use cr/ y/n as error on Ab, where cr is the standard 
deviation and n are the counts in each bin. 

4. Results and Discussion 

4.1. Oscillation amplitude VS magnetic field strength 

In Figure |2] we show the integrated amplitude A „ scatter plot 
(black dots), and Ab for each magnetic bin (red line), for the 
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Fig. 1. SDO-HMI mean magnetograms of (a) AR 1 1 166 (b) AR 
11316. See Table[T]for more details. 



ySARs shown in Figure [T] Hereafter, we refer to the As plot as 
the amplitude profile. The amplitude shows a decreasing trend 
up to a plateau, that is almost always reached in the range 0.5-1 
kG (as in the case of Figure I2J3). This feature suggests that the 
amplitude is independent of the magnetic field strength above a 
threshold value, which is different for each ySAR and each polar- 
ity in the same ySAR. 

For every polarity present in the 12 /3ARs analyzed, we com- 
puted the oscillation amplitude reduction as the ratio d = 
A1000M25, where A 1000 and A25 are, respectively, the mean os- 
cillation velocity around B - 1000 G, and that around B = 25 G, 
the minimum observable strength. Averaging the d values over 
all the ySARs, we obtained d - 0.54 + 0.06 , which agrees with 
the oscillation amplitude red uctions in magnetic environments 
quoted in the liter ature (e.g. Woods & Cramlll98ll iLites et aP 
119821: iBraun et al.| [l987: Tar bell et al.lll988h . 

The plateau in the amplitude profiles for high B values, 
is intriguing. By visual inspection, we verified that it takes 
place almost exclusively in the /3AR spot umbrae, as can be 
expected given the high magnetic field values. We can exclude 
an instrumental effect, since the magnetic strength threshold 
changes with each y6AR and even with the polarity in the same 
;6AR, and we are well above the 20 m/s error on the velocity 
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Table 1. Details of SDO-HMI data used. The fourth column 
shows the leading to trailing polarity area ratio. 



600 



/3AR number Observation Date Location 



S,/S, 



AR11166 


2011.03.08 


Nil WOl 


0.88 






711:00-14:00 




AR 11283 


2011.09.05 


N13 WOl 


0.86 






117:00-20:00 




AR 11302 


2011.09.28 


N13 W06 


0.96 






T17:30-20:30 




AR11319 


2011.10.15 


N09 E02 


0.95 






T07:00- 10:00 




AR11316 


2011.10.15 


S12 EOl 


0.90 






T07:00-10:00 




AR11330 


2011.10.27 


N09 E04 


0.83 






T18:30-21:30 




AR11338 


2011.11.06 


S 14 WOO 


0.92 






120:00-23:00 




AR11341 


2011.11.11 


N08 WOO 


1.05 






107:00-10:00 




AR 11362 


2011.12.03 


N08 W06 


1.10 






T20: 00-23: 00 




AR 11367 


2011.12.07 


S18 W05 


0.90 






T20:00-23:00 




AR 11375 


2011.12.13 


N08 W03 


0.82 






T19:00-22:00 




AR11389 


2012.01.03 


S22 W05 


0.97 






716:00-19:00 





measurements due to the HMI filter transmission profiles 
(Fleck et al. 2011). It seems that the oscillation amplitudes are 
not completely suppressed even in the strongest magnetic fields, 
and are not dependent on the ma gnetic field strength a nymore. 
Conversely to the simulations of Cattaneo et al.l ( l2003h . it may 
indicate a saturation effect of the kinetic energy in the strong 
field regime, and shed new light onto the mechanisms that 
co-operate in the reduction of acoustic power in the lower 
atmosphere. 



4.2. The leading-trailing asymmetry 

To investigate the relations between local oscillation amplitude 
and magnetic field strength, we compute the midpoints of the 
amplitude profiles. Specifically, we draw horizontal lines at con- 
stant Ab values, proceeding down in steps of 10 m/s. We com- 
pute the abscissa of each midpoint on the segment that intersects 
the amplitude profiles. If these profiles were symmetric, the mid- 
points should be on the B - vertical line. In order to compare 
the amplitude profiles from different ySARs, we normalize them 
to their average value in the B = +25 G bins. We remind 
that, for solar cycle 24, the leading polarities are negative in the 
N hemisphere and positive in the S hemisphere. Figure|3t shows 
the midpoint lines for all the 12 /?ARs considered in this study. 
Northern and Southern ySAR families appear well defined and 
separated. Moreover, as can be seen in Figure |3] the amount of 
asymmetry, defined as the distance of each midpoint from the 
vertical line through B - 0, increases for stronger fields. In the 
N hemisphere (in blue) it extends up to |B| ^ 230 G, while 
in S hemisphere (in red) it is less apparent, and only reaches 
|B| ^ 130 G. The amplitude profiles of Northern /3ARs are biased 
toward positive polarities, i.e. the trailing, while the Southern 
jSARs are biased toward negative polarities, i.e. the trailing again. 
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Fig. 2. In panel a) the integrated velocity oscillation amplitude 
for ARl 1166 is shown; the same for AR1136 in panel b). The red 
overplots correspond to the average in each magnetic field bin. 
In the inset of panel b) amplitude distribution in the 25 < B < 50 
G bin (black line) and gaussan fit (red line) are shown. 



To visualize differently the asymmetry, we can fix a value |B| 
for the magnetic strength, and compute the amplitude difference, 
i.e. the difference between the respective averaged amplitudes 
for both polarities in each /3AR, namely A|b| - A_|b|. In Figure 
[3J5 we show these differences versus the absolute value of mag- 
netic strength for all jSARs of Table [T] This plot confirms that 
the asymmetry is lower where the field is weaker At 25 G the 
difference is about 10 m/s and typical values for velocity am- 
plitudes up to 600 m/s; at 1000 G the difference rises to over 
50 m/s, while the amplitudes are ~ 200 m/s. Both in N and S 
hemispheres the oscillation amplitude in the trailing polarity is 
higher than that in the leading polarity. This fact implies that 
A|B| > A_|B| in the N hemisphere and viceversa in the S hemi- 
sphere, and that the amplitude reduction is less efficient in the 
trailing polarity. 

Such an asymmetry suggests that the amount of oscillation 
reduction within magnetic environments is not a function of 
the local magnetic field strength only. If this were the case, we 
would observe the same oscillation amplitudes for B and -B 
magnetic fields. This would imply both the plots of Figure [3] to 
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Fig. 3. a) Midpoints for normalized velocity oscillations ampli- 
tude, b) Difference between oscillation amplitude for both po- 
larities against the magnetic strength. Bars represent errors. Red 
line with diamonds, refers to Southern ySARs, blue line with plus 
signs to Northern ySARs. 



lie on the line through G (a) and m/s (b), respectively. 



The simulations of iFan et al l (fT99l about area asymmetries 
of ySARs, predicted a pronounced difference in the field strength 
between the leading and the trailing polarities, and ensuing 
differences in their fragmentation and lifetimes. This suggests 
that the asymmetry we found could derive from non-local 
properties, i.e. on the whole magnetic environment topology 
and even its surroundings. 

We therefore co mputed the area a symmetry in our magne- 
tograms following lYamamotol i20 1 2|) . with a threshold B > 500 
G. The results are reported in the fourth column of Table [1] We 
found that, on average, the area of the leading polarity is ~90% 
of the trailing polarity area, regardless the hemisphere which the 
jSAR belongs to. This quantifies how much a trailing polarity 
is spread with respect to its leading polarity. We can speculate 
that a relation may exist between the area and the oscillation 
suppression, due to a different interaction of acoustic waves 
with sparser or denser active region polarities. 



4.3. Oscillation amplitude VS magnetic field inclination 



The case ofAR11166 

To support the findings reported in sec. 14.2] we studied the de- 
pendence of the oscillation amplitude on the magnetic field in- 
clination angle in an ySAR. For this purpose, we considered 
the AR11166 (Figure la) as a case study. We used the LOS in- 
clination maps pr ovided by the VFIS V (Very Fast Inversion of 
the Stokes Vector, iBorrero et al.ll201 H) inversion of HMI vector 
magnetic field data. 

In Figure |4] we show the magnetic field inclination map of 
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Fig. 4. Inclination map of AR 1 1 166. The leading polarity (nega- 
tive) ranges from 90° to 180°, and the trailing polarity (positive) 
fromO° to 90°. 

AR11166. The inclinations retrieved by the VFISV span the 
range 10° < 6 < 172°. Following the most used sign convention, 
the leading polarity (negative) ranges from 90° to 180°, and the 
trailing polarity (positive) from 0° to 90°. As expected, the inner- 
most umbrae are mostly line of sight, as the /?ARs were selected 
as close as possible to the disk center 

VFISV also provide s the error erg associated to the inc lina- 
tion, for each pixel (iBorreroet al.ll201lt jPress et al.lll986l) . To 
avoid polarity flips due to the inversion error, we discarded all 
those pixels where 9 + 3crg causes the field to flip its inclina- 
tion from positive to negative (or vice-versa). We stress that 
VFISV is a Milne-Eddington code, therefore blind to any B 
and V line of sight gradients, which produce asymmetries in the 
Stokes profi les. This happens most often in regions filled by 
weak fields dViticchie et alJl20lll IViticchie & Sanchez Almeida 
l20Tlh . Moreover, the inversion of low signal to noise Stokes 
profiles is usually p roblematic (see iLandi DegITnnocentill992l : 
iBorrero et al.ll20I U for more details). For noisy Stokes profiles, 
and therefore low B values, VFISV (and an y inversion code) 
tends to be biased toward ~ 90° inclinations dBorrero & Kobell 
I20TTI) . For these reasons and to be consistent with the thresh- 
old used in |3] we rejected the pixels with B < 25 G and 
75° < < 105°. After this selection, about 17% of the initial 
pixels were selected for the analysis. 

We plotted Aw and its average in each 9 bin Ag =< A„. >g 
against the magnetic field inclination 9 (shown in Figure|5]l. For 
ease of comparison, in the upper panel of Figure |6] we show 
again Ag versus 6 in the range [10°, 75°]. Both Ag smoothly de- 
crease from ^ 250 m/s to ~ 140 m/s, from almost horizontal to 
almost vertical fields, respectively, but the oscillation amplitude 
reduction is more effective for the leading polarity (red curve) 
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Viticchie et al. 


|201l|). 



S 100 



100 

Inclinatbn |deg] 

Fig. 5. Integrated velocity oscillation amplitude within ARl 1166 
as a function of the magnetic field inclination 6. The thick red 
line corresponds to the average amplitude Ae -< A „ >g in each 
6 bin {A9 - 1°). The vertical dashed Unes mark the threshold 
6i<75°OR6»> 105°. 




30 40 50 
Inclination [deg] 

Fig. 6. Upper panel: Ag and Ag versus 6 in the range [10°, 75°]. 
The blue line represents the average amplitude of the positive 
(trailing) polarity: Ag. The red line represents the average am- 
plitude of the negative (leading) polarity after a remapping op- 
eration: 6 — > 180° - ft Ag. Lower panel: AA versus (see the 
text). Vertical bars represent the error, which were computed for 
each Ag with the same method described in section [3] and then 
summed up. 



than for the trailing polarity (blue curve). 

In the lower panel of Figure |6] we plot the oscillation amplitude 
asymmetry AA = - Ag VS 0. AA is constant (^ 15 m/s) 
within the errors up to ^ 45°, then AA drops to very few m/s. 
The inclination analysis shows that a ^ 15 m/s asymmetry exists 
between the polarities of ARl 1166. In particular, the velocity 
oscillation amplitude is enhanced in the traihng polarity with 
respect to the leading. 



and we recall that a VFISV hypothesis is that the magnetic field 
is constant within th e pixe l. 

ISchunker & Callvl (|2006|) and IStangahni et"al1 (1201 Ih have 
demonstrated that, due to mode conversion, there exists a pre- 
ferred angle at which the power is significant larger In our case, 
we indeed focus on the amplitude asymmetry between the two 
polarities, instead of considering the proper amplitude depen- 
dence that may be affected by biases due to the magnetic field 
strength. 

5. Conclusions 

In this work we focused on the velocity oscillation amplitude 
reduction in /JARs, reaching the following conclusions. 

The oscillation amplitude reduction found in magnetic envi- 
ronments of /?ARs is 0.54 + 0.06, which agrees with the values 
quoted in the literature (e.g. Woo ds & CramI I198U iLites et all 
ll982HBraun et al.lll987tlTarbell et al..,1988lK 



The five-minute amplitudes are not completely suppressed in 
the strongest magnetic fields of the spot innermost umbra. It 
would be very important to understand why a plateau appears in 
the oscillation amplitude profiles for strong fields. 

There exists a leading-trailing polarity asymmetry in /JARs. The 
asymmetry suggests that the reduction in oscillation amplitude 
does not depend on the field strength only, but may depend also 
on non-local conditions, such as the area on which the field 
spreads, for instance. 

The trailing polarity systematically shows a higher oscillation 
amplitude than the leading polarity, regardless the hemisphere 
which the ySAR belongs to. 

The plot of the velocity oscillation amplitude as a function of 
the magnetic field inclination confirms such an asymmetry in 
ySARs. Apparently, the asymmetry does not evidently depend on 
the inclination. 

We took advantage from HMI full-disk data at a 45 s time 
cadence, and at 1 arcsec spatial resolution. The analysis of these 
data revealed a possible saturation of the oscillation amplitude 
reduction for strong B, and an asymmetry in such a reduction for 
leading-trailing polarities. These results have to be accounted for 
in modeling the power reduction in magnetic environments, the 
emergence and the evolution of ySARs. 
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There is no evident relation between the asymmetry AA 
and the inclination 9. For any < 45° the oscillation am- 
phtude asymmetry is positive and around 15 m/s. We inter- 
pret the drop at ^ 45° as the combined effect of the VFISV 
preference to associate weak fields with larger inclinations 
(BoiTero & Kobel IQl l), and of the amplitude difference to 
be small for weak fields (see Figure ^p). Also, an incor- 
rectly retrieved weak B may result due to unresolved mag- 
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